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1. Introduction 

The renormalised QCD coupling a s = g 2 /(4n) is one of the fundamental parameters of the 
Standard Model of particle physics and can not only be determined from high-energy experiments 
but also estimated directly in lattice QCD simulations (see, e.g., for a recent review). 

The actual values of a s depend on both the renormalisation scheme (including the number of 
active flavours) and the scale considered. Given a renormalisation scheme S the dependence of af 
on the scale pL is controlled by the renormalisation group through 

,,2 d ccf(n 2 ) S( S , f a f\ l+2 

where j3 s is the beta function defined in that scheme. Its asymptotic expansion is known up to four 
loops in the MS scheme [Q] and up to three loops in various MOM schemes [^]. The first two co- 
efficients are renormalisation-scheme independent and can be looked up, e.g., in [Q] and references 
therein. The number of physically admissible renormalisation schemes (RSs) is unlimited, and so 
is the number of running couplings. The MS scheme, with the underlying use of dimensional reg- 
ularisation, is currently the most widely used RS in the analysis of high-energy experimental data. 
Such experiments are usually performed at different scales /I,, but through Eq. ( |l.l[ ) the different 
values of a^* s (jUf) are related to each other. 1 In order to keep track of different initial parametri- 
sations (ju,-, a^ s (/x,-)) a scale-invariant parameter A MS is introduced. This can be done for any RS. 
Once the A parameter is known in one RS, a one-loop calculation is sufficient to obtain it in any 
other scheme (see, e.g., [|7|]). 

Various nonperturbative studies in the past, in particular in lattice QCD, have provided esti- 



mates for a^ AS (Mz) or A MS for different numbers of flavour (see, e.g., [(8|-11]). There, different 
RSs, and thus different nonperturbative definitions of a s , have been used and the results roughly 
agree with what has been found in experiment. Our objective here is to show that, in future, good 
results for a s MS (Mz) or A MS may be expected using lattice QCD in Landau gauge. 

The definition of the running coupling we employ here is a nonperturbative one that has been 
first presented in the context of introducing a solvable systematic truncation scheme for the Dyson- 
Schwinger equations of Euclidean QCD in Landau gauge [|l2|]. This running coupling, which we 
call af 1 ™ in what follows, is defined in a MOM scheme and has its seeds in the ghost-gluon vertex 



in Landau gauge. To be specific, we use [ 12] 



ar M (q 2 ) := «r M (M 2 )Z(<7 2 ,MVV,M 2 ) d-2) 

which defines a nonperturbative running coupling that enters directly into the DSEs of QCD [|i~3|]. 
Z and / are the renormalised dressing functions of the gluon and ghost propagators, 

Df v { q \^) = 8 ab (8, v - q ^)^p- and G«V, M 2 ) = . (1.3) 



'in the literature it is common practice to evolve data to a^ s at the Z-Boson mass (see, e.g., [g, Fig. 23] for a nice 
illustration.) A recent compilation of a^(Mz) values has been given at the ICHEP conference resulting in the world 
average ce^Mz) = 0.1175 ±0.0011 [&. 
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In a lattice regularised theory, the corresponding bare dressing functions Zl and Jl are related to Z 
and J, renormalised at some sufficiently large q 2 = /I 2 , by 

Z(q 2 ,ii 2 )=Z- l {ii 2 ,a 2 )Z L {a 2 ,q 2 ) and J(q 2 ,H 2 ) =Z^\ i i 2 1 a 2 ) J L (a 2 ,q 2 ) (1.4) 

where Z3 and Z3 are the respective renormalisation constants and a is the lattice spacing. When 
considering I /a as our lattice UV cutoff, i.e., for a — > 0, then 2 

g 2 W Z 3 (/iV 2 )Z 2 (/i 2 ,a 2 ) o^O MOM/ 2x Q 5) 

47T Z^a 2 ) 1 ^ ^ 

With discretisation errors of 0(a 2 ) we can thus write 

«M° M (<? 2 ) = sM zdq 2^2 )J 2 {q 2y ) + 0{a 2^ (L6) 

where g 2 (a) is the bare coupling at the lattice cutoff scale 1 /a. It is this form of the running 
coupling which we will use below. 

2. Details of the numerical simulation 

The preliminary results described here were obtained on both zero and two-flavour SU (3) 
gauge field configurations. The quenched configurations were thermalised using the standard Wil- 
son gauge action at several values of j3 = 6/g 2 (a). We applied update cycles each consisting of 
one heatbath and four micro-canonical over-relaxation steps. The unquenched gauge field con- 
figurations were provided to us by the QCDSF collaboration. They used the same gauge action 
but supplemented it by Nf = 2 clover-improved Wilson fermions at various values of the hopping- 
parameter K. For details on the choice of /3- and fc- values we refer to Tab. [| All gauge con- 
figurations were fixed to Landau gauge with an iterative Fourier-accelerated gauge-fixing algo- 
rithm [18]. 3 For the stopping criterion we chose max x Tr [(V iU A^ iU )(V jU A v ^) t ] < 10~ 13 . The 
fields A x u =A fl (x + fl/2) are the lattice gluon fields given here in terms of gauge-fixed links U x ^ 
by the mid-point definition 

A^x + A/2) := jA-fe ~ U ln) -^r T < U ^- U ^ ■ 

This is accurate to order 0(a 2 ). On each such gauge-fixed configuration the momentum-space 
gluon and ghost propagators were measured. On the lattice, these are defined as the Monte Carlo 
averages 

Df v (k) = (A^(k)A b v (-k)) u and G ab (k) = I (v (M-^, e^A 



2 Note that, in Landau gauge, as shown long ago Jl4[] , the ghost-gluon vertex is regular and finite to any order in 
perturbation theory. Its renormalisation constant can thus be set to Z\ — 1 . Numerical evidence that this is also valid 
nonperturbatively has been provided in various investigations (see, e.g., jl5|-[l7|]). 

3 Note that for the range of momenta studied here the Gribov ambiguity is irrelevant as verified numerically in [n9l|. 
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p 


K 


latt. 


a [fm] 


#conf. 


P 


latt. 


a [fm] 


#conf. 


5.25 


0.13575 


24 3 x 48 


0.084 


60 










5.29 


0.13590 


24 3 x 48 


0.080 


55 


6.20 


32 4 


0.063 


30 


5.40 


0.13610 


24 3 x 48 


0.070 


62 


6.40 


32 4 


0.048 


50 


5.40 


0.13640 


32 3 x 64 


0.068 


57 


6.60 


32 4 


0.037 


60 


5.40 


0.13660 


32 3 x 64 


0.068 


30 


6.80 


32 4 


0.029 


46 



Table 1: Number of Nf = 2 (left) and Nf = (right) configurations used. To set the lattice spacings we use 
values for (ro/a) as provided by QCDSF || ^] and assume ro — 0.467 fm. 



where =A a ^T a are the Fourier-transformed gluon fields, M is the lattice Faddeev-Popov operator 
in Landau gauge, and k -x = J\, Ink^x^/L^. For a definition of M and details on its inversion we 
refer to [ |l7[ [TjJ and references therein. The corresponding bare dressing functions, Zi and Ji, 



are then calculated by assuming a tensor structure for the lattice propagators as given in Eq. (|1.3|), 
but with the continuum momenta substituted by p^ik) = (2/a) sin(7Ck^/L li ) where the integers 
k^ G {— L ll /2+ 1, . . . 4 At tree level this tensor structure is exact for the Wilson gauge action 

and the Faddeev-Popov operator we use (applying periodic boundary conditions). Given the data 
for Zl and Jl, we then estimate the coupling constant by the product 

oc L (p 2 ) = 8 ^Z L {p\a 2 )J 2 L {p\a 2 ) (2.1) 

where g 2 (a) = 6/(5. The corresponding error is obtained from a bootstrap analysis. Note that 
unlike other lattice investigations where data for Z and J are usually separately renormalised, no 
renormalisation is done here. 



3. Preliminary results 

In Fig. H (left) we show our present Nf = 2 data on Gfc as a function of momentum. One 
sees that the data for different values of the input parameters j3 and K form a reasonably smooth 
curve, indicating that discretisation effects are small for the input we employ. The only noticeable 
exceptions are the highest momentum values of each data set, where a slight deviation becomes 
evident with decreasing a. As expected, the quark-mass dependence is small. 

Discretisation effects are more evident for our present quenched data, as shown in Fig. |] 
(right). In particular for j8 < 6.4, which corresponds to even smaller lattice spacings than we 
use for Nf = 2, the data at larger momenta tends towards larger values as a decreases. Perhaps 
surprisingly, this effect is in the opposite direction than that observed in the unquenched data. 

Since we use a tree-level improved definition of momentum (see above) the data at larger 
momenta is expected to be less affected by discretization errors than with the naive definition 
ap^ = 2nk^/L^. Nevertheless, if the data is fit to perturbative QCD, e.g., as described below, the 
momenta p^ considered must satisfy A ~ 250 MeV \p^\ <C %ja . For the cylinder-cut data 
[21] used here, this translates into an upper bound of p 2 = p^p^ <C 4n 2 /a 2 ~ 300 GeV 2 with the 



lattice spacings for Nf = 2 in Fig. ^. This leaves us a considerable range of momenta to work with. 



4 Of course, for both propagators the case k = (ki,k2,k^,k4) =0 has to be excluded. 
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Figure 1: Nf = 2 (left) and Nf = (right) data for a L (p 2 ) at different j3 (and k). We set r = 0.467 fm to 
assign physical units to p 2 . Approximate values for the lattice spacing a are given in fm. 



We observe, however, that for both the quenched and the unquenched data, strong discretisation 
effects appear at a 2 p 2 > 14 (not shown), where the data points (for fixed parameters) start bending 
downwards as a 2 p 2 is further increased. We therefore restrict the present analysis to momenta 
a 2 p 2 < 14 where no such effect occurs. 

Given the data for different parameters, we consider, for each data set, a range of fitting win- 
dows within which the data is fit to a perturbative expansion of our running coupling. Since, from 
perturbative QCD, the running of gi^om is known up to three loops for the ghost-gluon vertex [9], 
we could use the truncated af 1 ™ beta-function at 3-loop order from Ref. [J3j] to describe the run- 
ning of our data. However, the coefficients c\ and ci in the expansion of q s MOM in terms of a^ s , 
i.e. 

a™ = af s { l + caf+c 2 [eg®] 2 + c 3 [af®] 3 + ...}, (3.1) 

which are known, as a function of Nf [pi, 5 are all positive and of 0(1) for the cases considered 
here, causing q s M to run more rapidly with scale than does q^ s . The result is that, as one 
lowers the scale, the truncated running at a given order becomes problematic at higher scales for 
the MOM coupling than it does for the MS coupling. In order to opt for the safest version of our 
analysis from the outset, we thus perform the running using the intermediate q MS running at four 
loops, with matching from af 1 ™ to a s MS at the start, and then re-matching back to q MOM at the 
end. 

In practice, this means that for our fits we scan over a fine-grid interval of a 4 MS values at an 
arbitrary reference scale pL, running each such value to all the momenta p 2 considered with the 
truncated 4-loop running 6 relevant to Nf. At each p 2 the 4-loop value of af ls (p 2 ) is then related to 
the corresponding MOM value q^ OM (/? 2 ) using the relation between the MOM and MS couplings 
given in Eq. (^J). For this we use the known coefficients c\ and C2 and set C3 = C4 = . . . = 0, 
an approximation which appears to become reliable for sufficiently large p 2 in our data sets. A 



5 For the ghost-gluon vertex those coefficients are roughly ci(0) ~ 4.23/jt and C2(0) ~ 36/7T 2 for the zero-flavour 
andci(2) ~ 3.67/7T and 02(2) ~26/n 2 for the two-flavour case 

6 To be specific, the running is performed using the exact analytic (implicit) solution of Eq. ( 1.1 ) corresponding to 
the four-loop truncated beta function. 
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(/3, k, a) = (5.40, 0.13660, 0.068) I— 

fit: c*4*-ioop — 6 — 




0.10 



20 40 60 80 100 120 

p 2 [GeV 2 ] 

Figure 2: Nf = 2 data for a L (p 2 ) at j3 = 5.4 and 
K" = 0.1366. The line represents the best fit to the 
data as described in the text. The fit window and the 
value of roA MS , as a result of that fit, are given too. 
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Figure 3: Fitted values for roA^ s as a function of 
bare mass am q for different /3 . We use K c values 
from [g fj(J] to set am q = (1/jc— l/)Q:)/2. 



^-minimisation of the scaled deviations of af lOM (p 2 ) from the a^p 1 ) data then determines the 
optimised af 1 (jU 2 ), from which the optimised A MS value is obtained using the conventional def- 
inition given, e.g., in Ref. We perform such a fit separately on the data for each (j8, K, Nf). 
An example of the resulting fit quality, together with the corresponding value for roA MS , is shown 
for the (j3, K,Nf) = (5.4,0.1366,2) set in Fig. ^[ Obviously, the data there is well described by the 
fit. The result shown there, roA = 0.60(1), is already in the ballpark of what is expected based 
on the results of, e.g., Refs. [^, 11]. In Fig. || all our present (as-yet-preliminary) values of roA MS 
for the Nf = 2 data sets are collected. Note that we have yet to fully investigate and quantify, and 
hence have not included in Fig. ^j, uncertainties due to (i) truncated running, (ii) the impact of pos- 
sible higher order, non-zero 03,04, . . . terms in the relation between the MS and MOM couplings, 
and (iii) the statistical uncertainty in the fit associated with that in the data. In particular, this will 
then allow us to extrapolate our roA MS values for the different parameters to the appropriate limits 
(a — > 0, K — *■ Kc) with realistic error estimates. 



4. Conclusions 

We have reported on first steps towards a determination of A MS in terms of lattice MC simu- 
lations of gluodynamics within the Landau gauge. Our method is based on the ghost-gluon vertex 
which in this particular gauge provides a nonperturbative running coupling in a MOM scheme de- 
fined solely in terms of the gluon and ghost dressing functions. Both these dressing functions can 
be calculated in terms of lattice MC simulations with good accuracy. 

Although our results are still preliminary, fits of our data to corresponding perturbative ex- 
pressions of our running coupling result in values of A^ s in the range expected. This suggests 
that a full estimate of the QCD parameter A within this framework is worth pursuing. Different 
systematic effects still have to be investigated before final conclusions can be drawn. In addition to 
those already noted above, lattice discretisation errors, in particular for Nf = 0, need further study. 
In the light of the results shown in Fig. ^, however, the approach appears quite promising. 
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